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NATIONAL  ADVTSOEY  OOMMTTTEE  FOR  AERONAUTICS 


TECBHICAL  NOTE  2222 

A  METHOD  FOR  THE  DETERMINATION  OF  THE  SPANWISE  LOAD  DISTRIBUTION 
OF  A  FLEXIBLE  SWEPT  WING  AT  SUBSONIC  SPEEDS 
B7  Richard  B,  Skoog  and  Harvey  H.  Brown 

SUMMARY 


A  method  is  presented  for  the  determination  of  the  apanwiae  load 
distribution  of  a  flexible  swept  wing  at  subsonic  speeds.  The  method 
is  based  on  a  relaxation  approach  utilizing  aerodynamic  loadings 
obtained  from  previously  published  work  (NACA  Rep.  921,  19k8)  baaed  on 
Weisainger’a  simplified  lifting— surface  theory  together  with  simple  beam 
theory.  The  solution  is  expressed  in  a  convenient  form  such  that  the 
amount  of  detailed  computing  involved  when  extensive  aeroelastic  cal¬ 
culations  for  many  flight  conditions  are  desired  is  reduced  to  that  for 
a  single  set  of  flight  conditions.  The  method  is  simplified  further  by 
an  abbreviated  solution  to  the  relaxation  process.  Sample  computing 
forBis  and  a  numerical  example  are  presented  to  illustrate  the  method. 


INTRODUCTION 


In  the  design  of  unswept  wings,  the  spanwise  distribution  of  aero¬ 
dynamic  load  usually  has  been  considered  to  be  Independent  of  structural 
deflections  since  torsional  deflecticjns  normally  are  negligible  at 
design  operating  speeds  (usually  well  below  the  wing  divergence  speed) 
and  bending  deflections  are  not  a  factor  for  zero  sweep.  On  a  swept 
wing,  however,  the  span  load  distribution  no  longer  may  be  considered  to 
be  Independent  of  structural  deflections  since  the  contribution  of 
bending  to  the  streamwise  change  in  section  angle  of  attack  can  become 
of  considerable  magnitude  as  the  sweep  angle  is  increased.  In  addition 
to  this  aeroelastic  effect  associated  with  wing  sweep  the  trend  toward 
higher  operating  dynamic  pressures  causes  factors,  which  were  previously 
negligible,  to  assume  Increased  significance. 

Several  methods  have  been  suggested  for  calculating  the  aerodynamic 
loading  of  the  flexible  swept  wing,  both  explicitly  (references  1  and  2) 
and  by  Iteration  (reference  3).  The  approach  in  each  of  these  methods 
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has  been  to  incorporate  aerodynamic  and  structural  theory  in  an  equation 
of  equilibrium  expressing  the  balance  existing  between  external  and 
internal  forces  on  the  wing.  To  simplify  the  mathematics  involved; 
matrix  notations  have  been  employed  in  two  of  these  methods  (references 
1  and  2).  Although  two  of  the  methods  (references  2  and  3)  indicate  how 
aerodynamic  theory  other  than  strip  theory  could  be  employed,  it  is^ 
apparent  that  such  application  increases  the  mathematical  difficulties 
a  great  deal  sc  that  the  simplest  aerodynamic  theory  has  been  used  in  all 
the  methods.  In  addition,  although  solution  of  the  problem  by  use  of 
matrices  has  the  advantage  of  simplicity  in  reducing  the  necessary  com¬ 
putations  tc  a  routine  form,  the  engineer  unversed  in  this  mathematical 
tool  encounters  a  loss  in  physical  appreciation. 

The  present  method  arose  from  the  effort  to  fill  the  need  for  a 
mathematically  si33:5)le  approach  which  can  yet  include  the  aerodynamic 
refinements  contained  in  lifting— line  or  lifting— surface  theories.  It 
was  desired,  also,  tc  separate  the  aeroelastic  effects  associated  with 
the  various  rigid-^ing  loadings  sc  that  it  would  not  be  necessary  tc 
perform  detailed  calculations  for  each  new  set  of  flight  conditions . 

The  method  of  this  report  is  based  on  relaxation  concepts  and  employs 
aerodynamic  span  load  distributions  from  previously  published  work 
(reference  4)  based  on  Weiss inger^s  simplified  lifting  surface  theory 
together  with  structural  deflections  found  from  simple  beam  theory. 

San^le  computing  forms  and  calculated  effects  for  an  example  wing  also 
are  presented. 

SYMBOLS 


A  wing  aspect  ratio 

a  distance  from  elastic  axis  to  section  quarter  chord  measured 

normal  to  elastic  axis  (positive  measured  forward),  inches 
(See  fig*  1*) 


av 


ratio  of  net  aerodynamic  force  along  the  airplane  Z  axis 

(positive  when  directed  upward)  to  the  weight  of  the  airplane 

wing  span  measured  normal  to  plane  of  symmetry,  inches 

section  chord  parallel  to  the  plane  of  symmetiy,  inches 

section  chord  normal  to  the  elastic  axis,  inches 

average  section  chord  parallel,  to  the  plane  of  symmetry.  Inches 

/c^  d7 


mean  aerodynamic  chord 


/  c  dj 


,  inches 


# 


u 


» 


section  pitching-inoment  coefficient  due  to  camber  as  usually- 
defined  f section  pitching  moment^ 


qc‘ 


section  pitching-moment  coefficient  due  to  camber  as  used  in 
this  report  (associated  with  camber^  chord  lengthy  and 
dynamic  pressure  normal  to  the  wing  quarter-chord  line) 


wing  lift  coefficient  ^ 

J 

rate  of  change  of  lift  coefficient  with  angle  of  attack  of 
root  section  at  plane  of  S3nimietry 

section  lift  coefficient 


section  lift  coefficient  from  additional-type  loading 

section  lift  coefficient  from  basic— type  loading  associated 
with  built-in  structural  twist 


total  change  in  section  lift  coefficient  due  to  structural 
deflection 


change  in  section  lift  coefficient  due  to  structural  deflec¬ 
tion  associated  with  additional  loading 

change  in  section  lift  coefficient  due  to  structural  deflec¬ 
tion  associated  with  basic  loading 

change  in  section  lift  coefficient  due  to  structural  deflec¬ 
tion  associated  with  torsional  moment  due  to  camber  loading 

change  in  section  lift  coefficient  due  to  structural  deflec¬ 
tion  associated  with  inertia  loading 

Young* s  modulus  of  elasticity^  pounds  per  square  inch 

total  change  in  section  angle  of  attack  due  to  structural 
deflection^  radians 

change  in  section  angle  of  attack  due  to  structural  deflection 
associated  with  additional  loading,  radians 

change  in  section  angle  of  attack  due  to  structural  deflection 
associated  with  basic  loading,  radians 

change  in  section  angle  of  attack  due  to  structural  deflection 
associated  with  torsional  moment  due  to  camber  loading,  radians 
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change  in  section  angle  of  attack  due  to  structural  deflection 
associated  with  inertia  loading,  radians 

huilt—in  structural  twist  of  tip  section  with  respect  to  root 
section,  radians 

change  in  section  angle  of  attack  due  to  structural  deflection 
produced  hy  rigid-wing  loadings,  radians 

change  in  section  angle  of  attack  due  to  structural  deflection 
produced  hy  the  aerodynamic  loading  introduced  hy  deflection, 
radians 

modulus  of  elasticity  in  shearing,  pounds  per  square  inch 

distance  from  elastic  axis  to  section  center  of  gravity 

measured  normal  to  elastic  axis  (positive  measured  forward), 
inches 

moment  of  inertia  in  bending,  inches  to  the  fourth  power 
polar  moment  of  inertia,  inches  to  the  fourth  power 
section  lift  loading,  pounds  per  inch 

section  lift  loading  plus  section  inertia  loading  for  chord 
sections  defined  by  c*,  pounds  per  inch 

bending  moment  at  any  position  along  the  swept  span  in  a 
plane  normal  to  assiomed  effective  root  (see  fig.  l), 
inch— pound 

dimensionless  span  coordinate  [Y/(b/2)],  fraction  of  semispan 

free— stream  dynamic  pressure ,  pounds  per  square  inch 

wing  area  (total),  square  inches 

semispan  measured  along  elastic  axis,  inches 

torsional  moment  at  any  position  along  the  swept  span  in  a 
plane  normal  to  the  elastic  axis,  inch— pound 

span  coordinate,  inches 
(See  fig.  1.) 

torsional  moment  representing  combined  contribution  of  section 
lift,  inertia,  and  moment  to  the  torsional  loading  about  the 
elastic  axis  at  any  spanwise  station,  inch-pounds  per  inch 
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V  angular  change  in  slope  of  the  elastic  axis  due  to  bending^ 

radians 

9  rotation  of  wing  sections  normal  to  the  elastic  axis  due  to 

torsion,  radians 

A  angle  of  sweep  of  the  quarter-chord  line,  degrees 

(See  fig.  1.) 

"'^a  angle  of  sweep  of  elastic  axis  (positive  for  sweepback),  degrees 
(See  fig.  1.) 

W  total  airplane  weight,  pounds 

w  section  structural  weight  based  on  chord  sections  defined  by 

c*,  pounds  per  inch 

\  wing  taper  ratio  (  ^P-  chord  | 

Vroot  chord/ 


THEOEY 

The  method  discussed  in  this  report  is  developed  to  apply  to  the 
general  case  of  a  flexible  wing  with  built-in  structural  twist,  camber, 
and  with  structural  inertia.  The  essential  feature  of  the  method  is 
the  application  of  relaxation^  procedures  to  the  physical  problem  of 
determining  the  aerodynamic  span  load  distribution  for  the  flexible 
wing.  In  formulating  the  theory  in  the  subsections  of  the  report  which 
follow,  the  unknown  aerodynamic  span  load  distribution  for  the  flexible 
wing  expressed  in  general  terms  is  applied  to  the  wing  together  with 
the  known  load  distribution  due  to  inertia.  From  the  bending  and  tor¬ 
sional  deflections  associated  with  this  loading,  the  rotation  (or  aero- 
elastic  twist)  of  wing  sections  parallel  to  the  plane  of  symmetry  then 
is  derived.  An  implicit  equation  thus  is  obtained  for  the  twist  dis¬ 
tribution  of  the  flexible  wing  corresponding  to  the  final  equilibrium 
position  of  the  wing  imder  the  combined  aerodsoxamic,  elastic,  and 
inertia  forces.  To  solve  the  equation,  relaxation  methods  are  applied, 
resulting  in  a  series-type  solution  for  the  loading  change  contributed 
by  structural  deflection.  As  a  final  step,  the  lengthy  series— typje 
solution  is  converted  to  an  abbreviated  solution  which  approximates  the 
series-type  solution  very  closely. 


In  reference  5^  relaxation  methods  are  shown  to  provide  a  powerful  tool 
in  solving  redundant  problems  of  structural  frameworks,  electrical 
circuits,  and  vibratory  systems.  In  the  present  report  the  relaxation 
approach  employed  is  based  on  the  principles  of  that  reference  rather 
than  the  exact  procedures  outlined  therein. 
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In  applying  relaxation  methods  to  the  present  problem^  the  wing  is 
assumed  to  he  fixed  in  position  as  the  rigid-^ing  loading  is  applied  to 
the  wing  and  then  the  wing  is  allowed  to  deflect  under  the  applied  load. 
The  wing  then  is  fixed  in  the  deflected  position,  while  the  loading 
corresponding  to  the  af ore-^nentioned  deflection  as  found  from  rigid-wing 
theory  is  applied  to  the  wing.  The  wing  then  is  allowed  to  deflect 
a^gain.  In  this  way  successive  deflections  can  he  found  which  are  depend¬ 
ent  on  the  loading  corresponding  to  the  previous  deflection.  In  the 
fomulation  of  the  theory,  this  procedure  is  used  only  to  develop  the 
series-type  solution.  Actually  in  practice  the  ahhreviated  solution 
previously  mentioned  is  used.  In  order  to  apply  the  method  most  expedi¬ 
tiously,  available  methods  for  the  quick  determination  of  aerodynamic 
span  load  distributions  for  wings  of  arbitrary  twist  and  arbitrary  plan 
foim  (e.g.,  reference  4)  should  be  used.  With  any  load  distribution  so 
determined  and  a  knowledge  of  where  to  apply  the  load,  it  is  then  a 
simple  matter  to  calculate  the  amount  of  structural  deflection  (bending 
or  torsion)  due  to  that  load  using  well-known  methods  of  solution. 

Before  discussing  the  theory  in  detail,  it  is  desirable  to  establish 
the  conventions  which  are  used  throughout  the  report  with  regard  to  span 
load  distribution.  In  accordance  with  the  convention  of  reference  6  the 
spanwise  distribution  of  lift  on  the  rigid  wing  is  considered  to  be  sep¬ 
arated  into  an  additional  lift  distribution  and  a  basic  lift  distribution 
tion,  the  additional  lift  being  proportional  to  wing  angle  of  attack  and 
the  basic  lift  being  dependent  only  on  built-in  structural  twist.  In 
this  report,  the  spanwise  lift  distribution  on  the  flexible  wing  is  con¬ 
sidered  to  be  separated  into  (l)  the  additional  and  basic  loadings  as 
defined  above  for  the  rigid  wing,  and  (2)  an  aeroelastic  loading  defined 
as  that  due  to  the  section  angle— of— attack  changes  produced  by  structural 
deflection. 

The  axes  referred  to  in  the  report  are  shown  in  figure  1.  The  Y 
axis  is  the  reference  axis  for  all  aerodynamic  span  load  distributions 
with  the  lift  assumed  to  act  at  the  quarter-chord  line.  The  wing  is 
assumed  to  have  an  effective  root  perpendicular  to  the  elastic  axis  and 
located  at  the  intersection  of  the  elastic  axis  and  the  plane  of  sym¬ 
metry. 

In  the  following  discussion,  the  basic  theory  is  developed  in  the 
subsections.  Loading  on  Flexible  Wing,  Aeroelastic  Twist  of  Flexible 
Wing,  and  Evaluation  of  Aeroelastic  Integrals.  The  details  of  this 
development  are  described  in  the  following  steps: 

1.  In  the  subsection.  Loading  on  Flexible  Wing,  expressions  are 
presented  for  the  running  load  and  running  torque  on  the  flexible  wing 
in  terms  of  the  component  loading  involved. 
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2.  In  the  subsection,  Aeroelastic  Twist  of  Flexible  Wing,  a  gen¬ 
eral  expression  is  developed,  for  the  twist  due  to  structural  deflection 
based  on  the  unknown  aerodynamic  loading^  for  the  flexible  wing  and  the 
effect  of  Inertia  as  presented  in  the  first  section.  It  is  also  shown 
how  the  general  equation  may  be  broken  down  to  lessen  the  amount  of  com¬ 
puting  involved  when  extensive  aeroelastic  calculations  are  desired, 

3.  In  the  subsection.  Evaluation  of  Aeroelastic  Integrals,  it  is 
shown  how  the  equation  resulting  from  application  of  the  relaxation 
method  can  be  abbreviated  to  provide  a  simple  and  direct  means  of  evalu¬ 
ating  the  implicit  twist  equations  presented  in  the  previous  subsection. 

Based  on  the  background  developed  in  these  three  sections,  an  equation 
is  presented  in  a  fourth  subsection.  Determination  of  Span  Loading  for 
Flexible  Wing,  showing  how  the  span  loading  for  the  flexible  wing  is 
determined  from  the  various  loadings  involved.  The  application  of  the 
basic  theory  to  determination  of  liftr-curve  slope,  aerodynamic  center, 
and  divergence  speed  for  the  flexible  wing  also  is  discussed  in  subse¬ 
quent  subsections. 


Loading  on  Flexible  Wing 


The  loads  which  will  produce  bending  of  a  flexible  wing  are  the 
rigid-H^lng  aerodynamic  lift  (additional  and  basic),  the  inertia  load 
normal  to  the  wing,  and  the  increments  in  aerodynamic  lift  produced  by 
aeroelastic  twist.  The  lifting  load  per  unit  span  based  on  streamwise 
wing  sections  can  be  expressed  as 


I 


(rigld-^ing  loading)  +  (loading  produced  by  wing  deflection) 
+  (c^^cq) 


(1) 


If  it  is  assumed  that  the  effect  of  taper  on  sweep  is  small  enough  that 
Ac/ 4  laay  be  taken  equal  to  A^a  and  if  the  effects  of  inertia  or  dead 
weight  also  are  included,  the  lifting  load  per  unit  span  along  the  wing 
panel  can  be  expressed  as 


^ In  setting  up  the  expressions  for  aeroelastic  twist,  it  is  convenient 
to  think  of  the  entire  flexible  wing  loading  as  being  unknown. 
Actually  all  components  are  known  except  the  component  Introduced  by 
structural  deflection. 
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(c 


2^cq+c^^cq)  cos 


Aqq-h^i  cq  cos 


Aea-^A^ 


(2) 


The  aerodynamic  and  inertia  loadings  normal  to  the  plane  of  the 
wing  will  produce  a  torsional  loading  about  the  elastic  axis  if  their 
lines  of  application  do  not  coincide  with  the  elastic  axis,  ^ If  a  wing 
has  camber  there  will  be  another  torsional  loading  due  to  this 

report,  the  contribution  of  camber  has  been  taken  such  as  to  affect  only 
torsional  deflections  of  the  wing  since  it  was  felt  that  values  of 
pitching-moment  coefficient  for  sections  normal  to  the  quarter-chord 
line  would  more  likely  be  available  from  two-dimensional  tests.  In  so 
doing,  however,  it  is  assumed  that  simple  sweep  concepts  can  be  applied 
to  finite— span  wings  without  serious  error.  The  torsional  moment  repre¬ 
senting  the  combined  contribution  of  section  lift,  inertia,  and  moment 
to  the  torsional  loading  about  the  elastic  axis^  at  any  spanwise  station 
thus  may  be  expressed  as 


t^  =  (rigid-^ing  torsional  loading)  + 

(torsional  loading  produced  by  wing  deflection) 
(cx^cq+c-^^-j^cq )a  cos  Aea  +  Cjj^^c^q  cos"^  Aqq^ 

(ci  cqa  cos  A© a) 
e 


>  (3) 


+ 


Aeroelastic  Twlat  of  Flexible  Wing 


In  this  section  of  the  report,  expressions  for  the  rotation  (or 
aeroelastic  twist)  of  streamwise  sections  of  the  flexible  wing  are 
derived  in  terms  of  simple  integrals  of  shear  and  moment  based  on 


®The  torsional  loading  about  the  elastic  axis  at  t]  =  Hi  (where  t)i  is 
an  arbitrary  value  of  t^)  will  be  a  summation  of  loadings  in  a  plane 
normal  to  the  elastic  axis.  Since  this  plane  will  intersect  the  span- 
wise  line  of  application  of  the  aerodynamic  or  inertia  loadings  at 
some  point  (where  t]2  ^  ili)  an  error  is  introduced,  which  was 

ignored  in  the  analysis. 
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elementary  beam  theory.  The  general  case  of  the  twist  due  to  the  load¬ 
ings  on  the  flexible  wing  given  by  equations  (2)  and  (3)  is  discussed 
first,  followed  by  discussion  of  a  convenient  way  of  breaking  down  the 
general  expression  for  aeroelastic  twist  to  facilitate  computation  when 
extensive  aeroelastic  calculations  for  a  large  number  of  flight  condi¬ 
tions  are  desired. 


Twist  due  to  loading  on  flexible  wing.-  The  bending  moment  at  any 
point  along  the  span  due  to  the  flexible-wing  loading  will  be"* 


Iw  dqdTj 


From  elementary  beam  theory 


and 


(6) 


(7) 


For  the  case  of  the  swept  wing,  "both  bending  along  the  elastic  axis 
and  twisting  about  the  elastic  axis  cause  changes  in  the  streanrwlse 
angle  of  attack.  The  change  in  section  angle  of  attack^  due  to  struc¬ 
tural  deflection  for  arbitrary  bending  and  twist  may  be  expressed  as 


^ Again  a  slight  error  is  introduced  (see  footnote  3)  because  the  various 
loadings  lie  along  different  axes.  To  be  correct  the  loadings  should 
be  referenced  with  respect  to  the  elastic  axis  rather  than  Y  axis. 
The  effect  of  drag  on  bending  moment  also  was  neglected  in  the  anal¬ 
ysis  and  normal  force  was  assumed  equal  to  lift  force. 

^This  involves  the  assumption  of  a  straight  elastic  axis  for  the  unde— 
fleeted  wing  and  ignores  the  effect  of  taper  in  regard  to  the  rotation 
of  wing  sections  due  to  torsion. 
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e 


—  V  ain  Aea  +  cp  cos  Aqq 


(8) 


In  the  case  of  swept— back  wings,  the  two  terms  of  equation  (8)  will  be 
of  opposite  sign.  For  swept-forward  wings,  the  terms  are  of  like  sign. 
Substituting  as  required  in  equation  (8),  we  have 


e  =  —  s  sin 


El 


dl]  +  3  cos 


=  —  s®  ain 


^1.0 

El 


di]  + 


> 


(9) 


S^  cos 


flo  ^ 

GJ 


dti 


J 


If  substitution  is  made  for  and  t^  from  equations  (2)  and  (3)  and 

if  the  terms  are  arranged  in  such  a  manner  as  to  show  clearly  the 
various  types  of  loading  which  contribute  to  the  atreamwise  twist  of  a 
flexible  wing,  the  expression  for  e  becomea:® 


© 

Equation  (lO)  consists  of  nine  terms.  The  word  description  given  oppo¬ 
site  each  term  on  page  11  merely  gives  a  physical  explanation  for  the 
existence  of  the  term  for  the  convenience  of  the  reader. 
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€  =  —  3  sin  A, 


COS  AeaCiridT] 


+  3^  COS  Aea/  — - - 


-  3=  3in  „ 


+  s^  cos  A, 


T)  cos  Aead.T) 


^  ^ea^n 

+  3^  cos  Aea/  — ^ - ; -  dt) 

Jq  GJ 


.  33  3in  a, 


p  ,  /"^/l.o'^hAzdT) 

-  32  COS  Aea/  - — -  dr) 

do  LrJ 


AeadTidn 


+  3  COS  Aea/  - - - 

Jo  GJ 


[1  ]  tvist  due  to  lend¬ 
ing  caused  by 
additional  load  of 
the  rigid  wing 

[2]  twist  due  to  tor¬ 
sion  caused  by 
additional  load  of 
the  rigid  wing 

[3]  twist  due  to  bend¬ 
ing  caused  by 
basic  load  of 
rigid  wing 

[4J  twist  due  to  tor¬ 
sion  caused  by 
basic  load  of 
rigid  wing 

[5]  twist  due  to  tor¬ 
sion  caused  by 
c?tin  of  rigid 
wing 

[6]  twist  due  to  bend— 
ing  caused  by 
inertia  loads 


[7]  twist  due  to  tor¬ 
sion  caused  by 
inertia  loads 


[8]  twist  due  to  bend¬ 
ing  caused  by 
aeroelastic 
loading 

[9]  twist  due  to  tor¬ 
sion  caused  by 
aeroelastic 
loading 


It  is  evident  that  the  first  seven  terms  of  the  preceding  expres¬ 
sion  give  the  twist  due  to  various  types  of  loading  of  the  rigid  wing. 

The  computation  of  this  portion  of  €  is  fairly  straightforward j  however. 
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the  last  two  terms  (terms  [8]  and  [9])  include  (which  is  a  function 

of  € )  so  that  an  explicit  solution  for  these  terms  is  not  immediately 
apparent . 

Twist  associated  with  individual  loadings,—  As  already  stated^,  the 
purpose  of  the  method  presented  in  this  report  is  to  determine  the  aero¬ 
dynamic  spanwise  load  distribution  for  a  flexible  wing*  If  this  load 
di-stribution  is  desired  for  only  a  single  set  of  flight  conditions 
(involving  a  specific  combination  of  q,  and  y  the  ovei^ll 

effect  of  flexibility  is  exi)re8sed  by  equation  (lO)*  However,  normally 
it  is  desired  to  examine  the  aeroelastlc  effects  of  a  wing  over  a  wide 
range  of  q,  and  A^;  and  for  various  combinations  of  these  factors. 

Therefore,  it  is  usually  more  convenient  to  separate  c^^  appearing  in 
terms  [8]  and  [9]  into  the  components  c,  ,  c,  ,  c.  .  ,  and  c, 

^ea  ^©b  ^©Cia^  ^^A^ 

associated  with  ^^b^  ^o^  respectively.  The  advantage 

of  solving  the  problem  in  this  way  lies  in  the  fact  that  only  one 
detailed  computation  of  the  components  of  c^^  need  be  made.  This 
simplification  arises  from  the  fact  that  (l)  ^2©^  proportional 

to  Cj^,  (2)  c^^  and  c^^  ^  do  not  depend  on  C^,  and  (3)  c^^^ 

is  proportional  to  A^  (or  to  Cj^  for  a  given  q  and  W/s),  It  is 
possible  also  to  perform  the  calculations  for  "unit  built-in  twist  and 
unit  camber  effect.  The  terms  of  equation  (lO)  therefore  can  be 
separated  into  several  components  as  follows: 


1.  Twist  for  additional  loading  only.-  The  aeroelastlc  twist 

associated  with  the  additional  loading,  consists  of  terms  [1] 

^d  [2]  of  equation  (lO)  plus  the  portions  of  terms  and  [9J 

contributed  by  .  Since  the  additional  loading  is  proportional 

®  a 

to  Ci,,  can  be  put  into  the  following  form: 


6cl  =  Cl 


^sinA 


T>  r f  ^  Aeadndii ' 

1 .0  1 .0  \CL°av/ 


ea 


El 


dTj  + 


^ea 


Cj  c 

r  *  i  _  /a 

^  I-  .0 


f  il 

Jo 


4) 


os  AeadTl 


GJ 


dTi  - 


qCg^s^sin  A 


ea 


/ 

Jo 


/Clg  Cv 

Tif'‘  (5:;^  jcosAeadTldn 

1 .0-^  1 .0  \^L^av/ _ _  + 

El 


cos  A 


ea/ 

do 


/.o 


GJ 


ea^dTi 

-S2 -  dTi 


(11) 
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In  this  equation  the  terms  within  the  brackets  are  based  on  unit  Cj . 
Once  these  terms  have  been  evaluated,  may  be  found  for  any  value' 
of  Cl  simply  by  multiplying  by  C^.  ^ 

2*  Twist  for  basic  loading  only.—  The  aeroelastic  twist  , 
associated  with  the  basic  loading,  consists  of  terms  [3]  and  [4]  of 
equation  (lO)  plus  the  portions  of  terms  [8]  and  [9]  contributed  by 
basic  loading  la  proportional  to  built-in  twist, 

can  be  put  in  the  following  form: 


qcavs®  sin  Aga 


r  ^  1  •o-'  1 .0  Y^b^Cav/ 

El 


cos  Aeg^dqdTi 


drj  4- 


1.0 


qcav^s2  cos  Aga.  J' 


T)  / 


a 

'av 


cos  Aga^dii 


GJ 


dT) 


(12) 


In  this  equation  the  terms  within  the  "brackets  are  based  on  unit 
structural  twist  of  the  tip  section  for  a  given  spanwise  distribu¬ 
tion  of  twist.  Once  these  terms  have  been  evaluated  for  this  dis¬ 
tribution,  can  be  found  for  any  amount  of  twist  having  the 

same  distribution  simply  by  multiplying  by 

3»  l^ist  for  camber  loading  only.—  The  aeroelastic  twist  €c*id^> 

associated  with  the  torsional  moment  due  to  camber  loading,  consists 
of  term  [5]  of  equation  (lO)  plus  the  portions  of  terms  [8]  and  [9] 
contributed  by  ^lecmQ*  torsion  due  to  camber  depends 

upon  the  amount  of  cainber  (or  the  value  of  Cmo)^  can  be  put 

in  the  following  form,  provided  the  caniber  is  constant  across  the 
span: 


n 
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'HIq  0 « 01 


qcav^s^  cos 


COS  A, 


dT]  — 


qc„„s3  sin  A, 


C  Q-rr 


COS  -AggdT|dT] 


1Cav^s2  cos  Aga 


S'" 

’I  1.0 


®  Ct» 


—  cos  AeadTi  . 
f - dTij  (13) 


In  this  equation  the  terms  vithln  the  brackets  are  based  on  a  imit 
camber  equivalent  to  c^iq  =  0.01.  Once  these  terms  have  been 
evaluated,  ^cmo  found  for  any  amount  of  this  type  camber 

by  multiplying  by  c^/O.Ol.  In  the  case  where  a  spanwise  varia¬ 
tion  in  cairiber  is  employed,  the  same  procedure  adopted  for 
can  be  used;  that  is,  calculate  the  terms  within  the  brackets  for 
a  given  camber  distribution  and  then  use  c'mQ  at  some  representa¬ 
tive  section  as  the  scale  factor. 

4.  Twist  for  inertia  loading  only.—  The  aeroelastic  twist  ^Az* 
associated  with  the  inertia  loading,  consists  of  terms  [6  ]  and  [7  ] 
of  equation  (lO)  plus  the  portions  of  terms  TS  ]  and  [9]  contributed 
'by  .  Since  the  inertia  or  dead-weight  loading  is  proportional 

®Az 

to  Av,  6Ai7  can  be  put  in  the  following  form: 


^Az  =  CavS®sinA, 


1  *0  1  *0 


qcavs®  sin  A^a 


1 .0  av  av 


dTj  — 


dT]  + 


qcav^s2  cos  Aga 


^®Az  Cav  , 
GJ 


: —  cos  Agadil  1, 
-  dq  (14) 
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In  this  equation  the  terms  ifithin  the  brackets  are  based  on  a 
structiiral  weight  at  Ig.  Once  these  terms  have  been  evaluated, 
e^„  may  be  found  for  any  load  factor  simply  by  multiplying  by  A^. 
It  shoiild  be  remembered  that  in  combining  the  effect  of  Av  given 
by  this  eqmtion  and  the  effect  of  Cx  given  by  equation  (11),  it 
is  necessary  to  adhere  to  the  following  relation,  in  order  to  retain 
any  physical  significance: 


Crq 


Evaluation  of  Aeroelastic  Integrals 


Swent-back  wings.—  As  has  already  been  discussed,  terms  [8]  and 
[9]  of  equation  (lO)  are  not  immediately  solvable j  however,  as  shown  in 
the  appendix,  the  solution  of  these  terms  can  be  expressed  as  a  power 
aeries  in  q  of  the  following  type: 

A€(ti)  =  fi(Tj)q^+f2(Ti)q®  +f3(Ti)q‘*+  .  .  .  (15) 

Equations  (lO)  through  (l4)  can  be  expressed,  therefore,  in  series  form 
as; 


e(Ti)  =  +fo(n)q+fi(Tl)q®+f‘2(ti)q®  +  f3(Ti)q‘*  +  .  .  .  fnCn)^”  (l6) 


where  the  coefficients  are  determined  by  the  particuleir  requirements  of 
the  equations  (lO)  through  (l4)  in  mind.  In  this  series  (equation  (l6)) 
the  valiies  of  successive  terms  are  of  opposite  sign.  If  equation  (l6) 
is  divided  by  fo(Tl)qj  we  have 


e(Ti) 

fo(n)q 


1  + 


fl(Tl)  ^ 

fo(T>)  **  ^  fo(n)  ^ 


+  .  .  » 


(17) 


As  has  been  shown  by  0.  K.  Smith  In  an  unpublished  Worthrup  report,"^ 
succeeding  terms  of  the  series  are  very  nearly  related  by  a  constant  of 
proportionality  so  that  equation  (I7)  can  be  written  as: 


e(ii) 

fo(n)'3 


1  —  kq  + 


(18) 


\  similar  approach  also  is  contained  in  outline  form  in  reference  7 
under  the  section  titled  Effect  of  Wing  G?wist. 
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where® 


f2(n) 

fx(Ti)  ' 


etc. 


As  also  shown  "by  Smith,  the  series  of  equation  (l8)  represents  an 
expansion  of  l/l+kq  so  that  eCq)  can  he  written  as: 

<«> 

In  these  equations,  foCq)  correspcsnds  to  terms  [1]  through  [7]  of 
equation  (lO)  (or  the  first  two  terms  of  equation  (ll),  (12),  or  (l4), 
or  the  first  term  of  equation  (I3))  with  q  set  equal  to  unity.  The 
function  fi(q)  corresponds  to  the  twist  produced  hy  the  aerodynamic 
loading  obtained  from  fo(il)  hy  the  method  of  reference  4.  With  fo(Tj) 
and  fi(T])  determined,  the  twist  dlstrlhutlon  of  the  fleiihle  wing 
e(Ti)  can  he  quickly  determined  for  any  value  of  q  hy  equation  (19). 
Then,  having  e(T|),  it  is  a  relatively  simple  matter  to  get  cj^(n)  hy 

the  method  of  reference  4,  which  is  baaed  on  the  Weissinger  simplified 
lifting-surface  theory  and  is  generalized  to  permit  determination  of 
load  distribution  for  a  wing  of  arbitrary  plan  form  and  arbitrary  con¬ 
tinuous  twist  distribution.  The  reference  can  be  used  to  provide 

^or  either  a  constant  lift  analysis  or  a  constant  angle-of— 
attack  analysis. 

It  should  be  noted  that  the  series  represented  by  equation  (l6) 
will  diverge  at  some  value  of  dynamic  pressure,  depending  on  the  struc¬ 
tural  rigidity.  An  outstanding  advantage  of  equation  (19)  (in  addition 
to  being  brief)  is  that  no  such  mathematical  difficulty  will  be  encotm- 
tered  so  that  the  aeroelastic  effect  at  any  dynamic  pressure  can  be 
calculated. 


Q  ■ 

In  practice,  it  is  usually  sufficiently  accurate  to  determine  k  as: 

k  =  _ 

Since  curves  of  the  functions  fo(Ti)  and  fiCrj)  will  generally  not  be 
of  exactly  the  same  shape,  the  value  of  k  will  vary  somewhat  across 
the  span.  In  the  case  of  the  example  wing,  the  best  approximation  was 
obtained  by  using  the  value  of  k  at  ti=1.0  since  at  that  station  the 
twists  given  by  successive  twist  distributions  (as  evaluated  using  the 
series— type  solution)  were  almost  exactly  proportional.  It  should  be 
noted  that  the  shape  of  a  given  e  curve  as  found  from  the  series-^ype 
solution  will  differ  slightly,  in  general,  from  the  shape  given  by 
fo(T)).  It  has  been  found,  however,  that  the  differences  in  curve  shape 
encountered  do  not  affect  C2^(t))  significantly. 


«r 
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Svept— f oarward  vlngs.—  The  method  can  he  applied  with  approximately 
the  same  accuracy  to  a  swept— forward  wing  hy  rewriting  equation  (19)  so 
that  a  minus  sign  appears  in  the  denominator  of  the  function  contain¬ 
ing  q . 

The  equation  then  becomes 


G  (n)  =  ^  ^oU) 

1  —  kq 

the  minus  sign  arising  from  the  fact  that  the  series  for  a  swept-f orward 
wing  (equivalent  to  equation  (l8)  for  a  swept— back  wing)  is 
l+kq+k^q^+  .  .  .  which  is  merely  an  expansion  of  l/l-icq. 


Determination  of  Aerodynamic  Span  Loading  for  Flexible  Wing 


The.  preceding  sections  of  this  report  have  laid  the  background  for 
determining  the  aerodynamic  span  load  distribution  for  a  flexible  wing. 
From  equation  (l) 


I 

q 


+  Cj^C  +  c 


(21) 


From  equation  (19)  it  is  evident  that  the  distrihution  of  across 

the  span  can  be  written  as 

cigC(Ti)  =  f  [  fo(Tl)]  (22) 

For  q=1.0,  this  load  distribution  can  be  written  as 


[C7  ®U)] 


q=  1.0 


1+k 


f  [fn(Tl)] 


(23) 


If  solution  of  f  [fo(Ti)]  from  equation  (23)  is  substituted  in  equa¬ 
tion  (22)  and  the  resulting  expression  for  c^^c  is  then  substituted 
in  equation  (21),  equation  (21)  can  then  be  written 


I 

q 


q(l+k) 


.0 


(24) 


If  the  following  relation 


cr  =  cr  +01  +  Cl  , 

'■e  ‘■ea  *^6^  ‘•ec’ 


+  Cl 


“o 


®A, 


is  substituted  in  equation  (24)  and  if  the  various  terms  are  written  as 
loading  coefficients,  equation  (24)  becomes 
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(25) 

With  the  various  aeroelastic  loadings  evaluated  for  q=1.0,  it  is  a  simple 
matter  to  calculate  (^/qcav)(Ti)  for  the  flexible  wing  for  any  combination 
of  q,  and  kz* 

Stability  Parameters  for  Flexible  Wing 


Lift-curve  slope.—  The  lift-curve  slope  for  the  flexible  wing  can 
be  found  by  either  of  two  methods ,  depending  on  whether  the  viewpoint 
adopted  originally  was  that  of  constant  lift  coefficient  or  constant 
angle  of  attack.  Both  methods  are  presented  here  for  convenience: 

1.  Constant  lift~coefficient  analysis.—  The  aeroelastic  span 
load  distribution  resulting  from  a  constant  lift  analysis  is  a 
basic— type  loading^  which  yields  zero  lift  when  integrated.  In 
solving  for  the  aeroelastic  loading  by  the  method  of  reference  ^4-, 
the  angle  of  attack  of  the  wing  root  required  to  obtain  zero 

over-all  lift  also  is  obtained.  This  angle  of  attack  represents 
the  angle  through  which  the  wing  root  must  be  turned  in  order  to 
maintain  a  given  wing  lift  coefficient  for  the  flexible  wing  at  all 
dynamic  pressures.  If  the  lift— curve  slope  for  the  rigid  wing  is 
known,  the  lift— curve  slope  for  the  flexible  wing  can  be  found 
graphically  by  simply  laying  off  the  value  of  cx^.  in  the  proper 
direction  from  the  rigid-^ing  lift  curve  ( increasing  wing  angle  of 
attack  for  swept-back  wings)  at  the  Cl  for  which  the  value  of 
was  obtained.  This  procedure  is  indicated  in  the  following  sketch 
for  swept— back  wings: 
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The  value  of  cXp  is  made  up  of  (l)  the  contribution  associ-- 
ated  with  rigid-^ing  additional  loading  and  (2)  the  contribution 
associated  with  inertia  (or  dead  weight)  loading^  since  only  the 
contributions  associated  with  these  loadings  are  proportional  to  C^. 
These  contributions  are  given  by  the  values  of  cx^  associated  with 
(t|)  and  c?-  (il)>  respectively.  These  loadings  are  found  from 

equations  (ll)  and  (l4)  previously  presented.  Lift— curve  slope  for 
rigid  wings  is  given  directly  in  reference  k  for  a  wide  range  of 
plan  forms. 

2  .  Constant  angle— of --attack  analysis.-^  The  aeroelastic  span 
load  distribution  resulting  from  a  constant  angle— of-attack  analysis 
is  similar  to  an  additional— type  loading  and  yields  a  lift  when 
integrated.  For  a  swept— back  wing,  the  change  in  lift  resulting  from 
integration  of  the  aeroelastic  loading  associated  with  the  twist  of 
equation  (19)  can  be  written  as 

ACi,  =  ^  <^1,  (26) 

where 

lift  coefficient  resulting  from  integration  of  the  aeroelastic 
loading  corresponding  to  e(T]) 

lift  coefficient  resulting  from  integration  of  the  aeroelastic 
loading  corresponding  to  fo(T]) 
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The  lift  coefficient  for  the  flexible  wing  (swept-hach)  therefore 


Cl  CLq  ^Li 


where 


CLq  lift  coefficient  of  the  rigid  wing 

If  equation  (27)  is  divided  by  the  arbitrary  angle  cf  attack,  the 
lift— ciurve  slope  for  the  flexible  wing  can  be  written  as 


q 

1+kq  CLq 


where  the  subscripts  F  and  E  refer  to  the  flexible  and  rigid 
wing,  respectively.  In  the  case  of  a  swept-forward  wing,  equation 
(28)  becomes 


1 

l-toi  CLq 


Aerod.yTin.Tni  c  ceirber The  accepted  definition  of  aerodynamic  center 
of  a  rigid  ving  is  the  centroid  of  all  the  additional  loads  •  It  is  evi~ 
dent  that  on  a  flexible  wing  the  lift  increment  due  to  angle  of  attack 
includes  not  only  an  additional— type  load  (in  the  rigid-^ing  sense)  but 
also  a  varying  amount  of  aeroelastic  lift.  The  effective  aerodynamic 
center  of  a  flexible  wing  will  thus  differ  from  that  of  the  rigid  wing 
but  will  still  be  the  centroid  of  all  the  additional  loads.  In  stability 
analyses^  it  usually  is  customary  to  neglect  the  vertical  location  of  the 
centroid  since  the  effect  of  drag  on  stability  normally  is  negligible. 

The  varying  amount  of  aeroelastic  lift  is  made  up  of  (l)  the  camr- 
ponent  associated  with  the  additional  loading  and  (2)  the  conrponent 
associated  with  the  Inertia  (or  dead  weight)  loading^  since  only  the 
conrponents  associated  with  these  loadings  are  proportional  to 

These  components^  c?  (tj)  and  ci  (tj)^  respectively,  are  defined  by 

©a 

the  twist  distributions  given  by  equations  (ll)  and  (l^)  previously  pre¬ 
sented.  With  (t))  and  c^^  (t))  determined,  the  aerodynamic  cen- 

ter  of  the  flexible  wing  can  be  found  by  any  method  for  determining 
centroids,  remembering  that  the  chordwise  load  is  assumed  to  lie  along 
the  quarter— chord  line  of  the  wing.  The  same  method  is  applicable,  in 
general,  whether  the  analysis  adopted  is  for  constant  lift  or  constant 
angle  of  attack. 


BACA.  TN  2222 


21 


Divergence  Speed  for  Flexible  Wing 


The  ease  with  which  divergence  speed  can  be  obtained  from  ecjuations 
(19)  and  (20)  is  worthy  of  brief  mention  also.  If  equation  (I9)  is  dif¬ 
ferentiated  with  respect  to  q,  the  derivative  becomes 

^  ■  (I^  'o*"’  '3°) 

The  condition  for  divergence  is  that  the  derivative  (d/dq)  e(Tj)  must 
eqml  infinity.  To  satisfy  this  condition, 

1  +  kqj)  =  0  (31) 

where  q^  is  the  divergence  dynamic  pressure,  so  that 

<3D  =  -  I  (32) 


In  the  case  of  a  svept— forward  wing,  differentiating  equation  (20)  and 
proceeding  as  before  yields 


Id  =  I  (33) 

« 

Equations  (32)  and  (33)  show  the  familiar  fact  that  a  sw'ept—forward  wing 
will  diverge  at  some  positive  value  of  q  while  a  swept-^ack  wing  will 
1  not  diverge  at  any  positive  value  of  q. 


APPLICATION 


Computing  Eorms 


Computing  forms  to  aid  in  solving  for  the  change  In  span  load  dis¬ 
tribution  due  to  wing  flexibility  are  presented  in  table  I.  The  basic 
data  needed  to  perform  the  calculations  are  shown  in  table  I  (a).  Tables 
l(b),  l(c),  l(d),  and  l(e)  are  for  the  purpose  of  computing  the  functions 
(fQ(Ti)  and  fi(T]))  for  each  of  the  rigid-wing  loadings,  and  require  essen¬ 
tially  the  same  operations.  Any  differences  noted  are  merely  for  ease  in 
handling  the  computations  for  a  given  loading.  Table  l(f)  is  essentially 
the  same  as  table  l(c)  which  Is  for  basic— type  loadings.  The  only  dis¬ 
tinction  between  the  two  forms  is  that  table  l(f)  is  for  the  aeroelastic 
loading  introduced  by  wing  deflection  and  table  l(c)  is  for  the  basic 
loading  due  to  built-in  twist.  With  one  exception  these  forms  consist  of 
separate  calculations  of  bending  deflection  (columns  4  throu^  7)  an^  of 
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torsional  deflection  (columns  9  through  11).  The  single  exception  Is  , 
table  l(e)  for  camber  loading  (cVq^  which  does  not  require  bending  com¬ 
putations  for  obvious  reasons.  The  last  column  (column  12)  in  each 
case  is  the  sum  of  the  bending  and  torsional  deflections  in  terms  of 
streamwlse  angle— of— attack  change.  This  coliimn  yields  fo('n)  in  the 
case  of  tables  l(b),  l(c),  l(d),  and  l(e)  and  yields  fi(T|)  in  the 
case  of  table  l(f).  The  column  headings  are  either  self-explanatory  or 
are  explained  in  the  computing  instructions  following  the  tables.  In 
the  computing  instructions  operation  A  Integrates  the  running  load 
nonaal  to  the  wing  to  obtain  the  shear  at  designated  spanwise  stations. 
Operation  B  integrates  the  shear  curve  so  obtained  to  determine  the 
bending  moment  at  the  same  spanwise  stations.  Operation  D  integrates 
the  running  torsional  load  to  obtain  the  torsional  moment  at  the  same 
spanwise  stations.  The  summations  are  performed  as  indicated  from 
T)  =  1.0  to  T)  =  0.  Operation  C  integrates  the  m/EE  curve  to  obtain 
the  slope  of  the  elastic  axis  at  the  chosen  spanwise  stations.  Opera¬ 
tion  E  integrates  the  T/GJ  curve  to  obtain  the  twist  about  the  elastic 
axis  at  the  same  spanwise  stations.  The  summations  in  these  two  opera¬ 
tions  are  performed  as  indicated  from  q  =  0  to  t]  =  1.0.  The 
integrating  operations  have  been  set  up  In  accordance  with  the  trapezoi¬ 
dal  rule  for  approximate  integration.  The  spanwise  stations  used  in 
the  computations,  therefore,  should  be  of  sufficient  number  and  of 
adequate  spacing  so  that  the  integrations  will  not  be  subject  to  signif¬ 
icant  error. 


Nmerical  Example 


The  solution  procedure  indicated  in  the  preceding  section  of  this 
report  has  been  applied  to  a  relatively  flexible  example  wing  of  moder 
ate  sweep  and  high  aspect  ratio.  Compressibility  considerations  were 
neglected  in  the  present  example,  since  the  modifying  influence  of  com¬ 
pressibility  is  small  compared  to  the  isolated,  primary  influence  of 
dynamic  pressure.  The  geometry  of  the  example  wing  is  shown  in  figure  2 
together  with  curves  of  the  structural— stiffness  data  used  in  the  cal¬ 
culations.  The  wing  has  an  aspect  ratio  of  9.^3^^a  taper  ratio  of  0.42, 
and  a  sweep  angle  of  the  quarter— chord  line  of  35  •  indicated  in 
the  figure,  the  elastic  axis  is  located  at  38-percent  chord. 

The  rigid-wing  loading  curves  used  in  the  calculations  are  shown  in 
figure  3.  The  aerodynamic  loadings  (additional  and  basic)  were  obtained 
from  reference  4,  neglecting  compressibility.®  The  additional— tjrpe 


^Eor  the  reader  interested  in  including  conpressibility,  it  should  be 
noted  that  the  method  of  reference  4  accounts  for  compressibility 
effects  on  the  basis  of  linearized  compressible  flow  theory  so  that 
the  modifying  influence  of  compressibility  can  easily  be  included 
using  that  reference  by  merely  following  the  procedure  outlined 
therein. 
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loading  corresponds  to  an  over-all  wing  lift  coefficient  of  1.0  and  the 
basic  loading  corresponds  to  a  linear  twist  distribution  having  1.0° 
washout  at  the  wing  tip.  The  dead-weight  shear  distribution  was  obtained 
from  weight  estimates  for  the  wing.  The  abrupt  shear  changes  shown  are 
due  to  the  concentrated  masses  of  the  engines.  The  camber  loading  was 
obtained  assuming  constant  camber  across  the  span  equivalent  to 

Cjijq  =  —  0.01. 

Loading.—  Computations  for  the  example  wing  assuming  only  the  addi¬ 
tional  loading  to  exist  and  based  on  a  constant  lift  analysis  are  pre¬ 
sented  in  table  II  to  illustrate  une  of  the  computing  forms.  Table 
II (a)  presents  the  geometric  and  structural  parameters  for  the  example 
wing  for  specified  stations  along  the  semispan.  Table  Il(b)  presents 
the  calculations  for  the  angle-of -attack  redistribution  from  which  the 
function  fp(T|)  is  found.  The  function  fiCq)  is  found  from  the 
angle-of -attack  redistribution  calculated  in  table  Il(c).  In  table 
II (b),  the  loading  used  Is  the  additional  loading.  In  table  Il(c),  the 
loading  used  is  the  basic  loading  found  from  the  twist  distribution  cal¬ 
culated  in  table  Il(b),  The  functions  fo(^)  and  fiCri)  and  the  ratio 
are  plotted  in  figure  4  against  spanwise  station  tj.  The 
ratio  fi(Ti)/fo(n)  la  plotted  for  the  purpose  of  illustrating  the 
differences  in  cu2rve  shape  between  the  functions  fo(Ti)  and  fi(Ti). 

As  can  be  seen  from  the  figure,  the  distributions  of  twist  (as  defined 
by  the  shape  of  the  curves)  for  fo(Tl)  and  fi(Ti)  are  somewhat  dif¬ 
ferent  due  to  the  fact  that  fo(il)  was  determined  from  an  additional 
loading  and  fiCi])  was  determined  from  a  basic  loading.  In  spite  of 
the  difference  shown,  however,  the  assumption  of  proportionality  between 
successive  terms  of  the  series  is  sufficiently  close  to  provide  a  good 
approximation  since  basic- type  loading  is  affected  very  little  by  dif¬ 
ferences  in  curve  shape  such  as  shown  in  the  figxjre.  As  can  be  seen 
from  the  figure,  a  large  variation  in  fi(T))/fQ(Ti)  across  the  semispan 
can  exist  so  that  it  is  important  to  choose  the  value  of  this  ratio  at 
the  proper  value  of  As  stated  earlier,  the  best  approximation  to 

the  series— type  solution  was  obtained  for  the  example  wing  by  choosing 
the  value  of  fi(ii)/fo(il)  at  the  tip  station.  These  remarks  also  apply 
if  a  constant  angle— of -attack  analysis  is  adopted. 

The  aeroelastlc  loadings  associated  with  each  of  the  rigid-wing 
loadings  are  presented  in  figure  5  for  several  values  of  dynamic  pres¬ 
sure  as  obtained  from  the  constant- lift  analysis.  With  the  solution  in 
the  form  shown  in  the  flgiire,  the  computation  of  the  spanwise  load 
distribution  for  the  flexible  wing  can  be  found  relatively  simple  for 
any  set  of  flight  conditions  and  any  set  of  camber  and  twist  distribu¬ 
tion  similar  to  that  aas\imed  merely  by  combining  the  loading  for  the 
rigid  wing  with  the  proper  combination  of  basic  loadings  due  to  deflec¬ 
tion  as  indicated  by  equation  (25) .  As  has  already  been  pointed  out, 
such  a  computing  shortcut  is  made  possible  by  the  linear  aerodynamic 
and  structural  theory  of  the  analysis  which  renders  the  deflection 
loadings  (at  a  given  dynamic  pressure)  proportional  to  either  Cj^,  A^, 
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or  Cjn^y  depending  on  the  rigid-^ing  loading  involved.  The  aero- 
elastic  loadings  for  a  constant  angle-of -attack  analysis  are  not  pre¬ 
sented  since  the  proper  combination  of  aeroelastic  loadings  is  believed 
to  be  indicated  sufficiently  by  figure  5. 

Stability  parameters.—  Since  no  experimental  data  vere  available 
with  which  theoretical  span  load  distributions  of  a  flexible  wing  could 
be  compared^  the  validity  of  the  method  was  checked  to  some  extent  by 
comparing  predictions  of  lift— curve— slope  change  and  aerodynamic— center 
shift  with  those  predicted  by  the  method  of  reference  1. 

The  variation  in  lift— curve  slope  with  dynamic  pressure  for  the 
example  wing^  as  calculated  by  the  procedures  outlined  earlier  in  the 
report^  is  presented  in  figure  6^  neglecting  the  modifying  influence  of 
compressibility.  In  the  calculations^  a  rigid-^ing  lift-curve  slope  of 
0.071  obtained  from  reference  k  was  used.  In  the  figure^  a  curve 
showing  the  slope  variation  according  to  the  method  of  reference  1  also 
is  presented  for  comparison.  As  can  be  seen,  agreement  between  the  two 
methods  is  very  good. 

The  spanwise  shift  in  the  centroid  of  the  additional  loads  with  . 
dynamic  pressure  for  one  panel  of  the  example  wing  is  presented  in 
figizre  7  together  with  the  corresponding  shift  in  aerodynamic  center 
parallel  to  the  plane  of  symmetry.  At  a  dynamic  pressure  of  500  pounds 
per  square  foot,  the  spanwise  shift  is  shown  to  be  inboard  by  about 
6  percent  of  the  semispan.  The  corresponding  chordwise  shift  is  shown 
to  be  forward  by  about  20  percent  of  the  mean  aerodynamic  chord.  Good 
agreement  with  the  prediction  of  reference  1  is  shown  again. 


Ames'  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  July  31>  1950. 
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APPENDIX 

DERIVATION  OF  WING  TWIST  POWER  SERIES 


In  the  text  of  this  report,  the  streaniwlse  twist  of  a  flexible, 
swept  wing  la  given  by  equation  (lO)  as  a  summation  of  terms  due  to  the 
various  rigid-wing  loadings  and  the  aeroelastic  loading  introduced  by 
flexibility.  Equation  (lO)  can  be  summarized  as  follows: 

€(ti)  =  (Al) 


where 

€q(t])  the  twist  of  the  flexible  wing  due  to  the  total  rigid-wing 

loading  (as  given  hy  terms  [  l]  through  [7]  of  equation  (lO)) 

and 

A€(t))  the  twist  of  the  flexible  wing  due  to  the  aeroelastic  loading  (as 
given  by  terms  [8]  and  [9]  of  equation  (lO)) 

As  stated  in  the  text,  the  value  of  A€(ti)  in  the  above  equation 
cannot  be  evaltxated  directly  in  any  explicit  manner  so  that  a  method  of 
successive  approximations  must  be  employed. 

One  method  for  evaluating  €(t|)  is  to  adopt  the  obvious  iterative 
approach  and  compute  ^(tj)  by  successive  approximations  to  Ag(t|)  until 
sufficient  accuracy  is  obtained.  A  more  convenient  method,  however,  is 
to  apply  a  relaxation  procedure  to  the  determination  of  A^(t)).  In  this 
method,  the  wing  is  assumed  to  be  fixed  in  position  as  the  rigid-wing 
loading  is  applied,  and  then  the  wing  is  allowed  to  deflect  under  the 
applied  load,  resulting  in  the  streamwise  twist  distribution 
equation  (A1)  .  The  wing  then  is  fixed  in  position  again  and  the  rigid- 
wing  loading  is  removed.  The  aerodynamic  loading  corresponding  to 
Gq(t])  then  is  applied  to  the  wing  and  again  the  wing  is  allowed  to 
deflect  (in  accordance  with  the  applied  load  only),  resulting  in  a  new 
twist  distribution.  In  this  way  successive  twist  distributions  can  be 
found  which  are  dependent  upon  the  aerodynamic  loading  corresponding  to 
the  previous  twist  distribution.  The  final  twist  distribution  for  the 
flexible  wing  can  therefore  be  expressed  as 


e(Ti)  =  eQ(q)+Aei(T])+Ae2(Ti)  .  .  .  +  Ae^_^(Ti) 


Comparison  of  equation  (A2)  with  equation 
Ae(n)  =  Aei(Ti)+Ae2(n)  • 


(A1)  shows  that 

.  .  +  A^ii) 


(A2) 

(A3) 
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Now  since 


Aei(ii)  =  q  F  Cj  (tj) 


where 


C,  (t])  =  f  6o(ll) 
^ei 


€o(’l)  =  <1  ^ 


^0(1)  =  <1  fo(l) 


It  is  apparent  that  equation  (a4)  can  he  written  as 


Aei(n)  =  qL^fa(n) 


Similarly^ 


Ae2(T))  =  q  F  c,  (t)) 


where 


so  that^  with  equation  (A5)^  it  is  apparent  that 


AG2(t])  =  q  f2(T^) 


In  like  manner 


A€^(t,)  =  fjri) 


In  these  expressions^  the  notations  f 


and  F 


merely 


indicate  the  general  dependence  of  loading  and  twist  on  the  associated 


twist  and  loading,  respectively.  The  notation  fj;^ 
a  specific  function. 


indicates 
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It  should  he  noted  that  for  a  swept— hack  wing  each  new  twist  dis¬ 
tribution  calculation  hy  this  procedure  will  he  of  opposite  sign  to  the 
twist  immediately  preceding,  since  successive  aerodynamic  loadings  will 
exert  bending  moments  of  opposite  sense  due  to  the  bending-twist  rela¬ 
tionship  for  a  Bwept-hack  wing.  Slimming  the  successive  contributions 
to  Ae(Ti)  of  equation  (Al)  as  given  by  equations  (A5)^  (AT)^  and  (A8), 
we  have 


Ae(ii)  =  q2fi(Ti)  +  q®f2(Tj)  +  .  .  .  q^'^’^f 3a( n )  (A9) 


Substituting  in  equation  (Al)  for  eo(T))  by  the  relation  given  under 
equation  (A^f)  and  for  A€(tj)  as  given  by  equation  (A9),  it  can  be  seen 
that  the  twist  distribution  for  the  flexible  wing  can  be  expressed  as  a 
power  series  in  q  as  follows: 


+  q®f2(Tl)  +  .  .  .  q“ftt-i(Ti)  (AlO) 


% 
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Spanwise  station,  ^ 

Figure  2-  Sketch  of  example  wing  with  pertinent 
dimensions,  plan- form  parameters,  and  elastic 
properties. 


Dead-  weight  shear,  lb  Basic  load  coefficient  Additional  load  coef- 
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(a)  Additional. 


(b)  Basic. 


Spanwise  station,  y 


(d)  Pitching  moment  due  to  camber. 


Figure  3.-  Rigid-wing  loadings  for  the  example  wing. 


Figure  4.  -  Th 
the  ratio 


fated  fa, 


Aeroelastic  deflection  per  unit  Aeroeiastic  loading  coefficient  per  unit 

lift  coefficient,  (e^),  rad  aft  coefficient,  c,  c/c,  c„. 
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{Rigid  wing 
additional 

I 

q,  psi 


Spanwise  station,  7^ 

(a)  Additional  loading. 

Figure  5.  -  Aeroelastic  loadings  and  deflections  associ¬ 
ated  with  the  various  rigid -wing  loadings  for  the 
example  wing  as  calculated  for  several  values  of 
dynamic  pressure  using  a  constant- lift  analysis. 


NACA  TN  2222 


q,  pS! 
0  ^ 


1^,- 

§  o-* 

-Js 

I  5 

^  Cl 


/f^/. 


-.02 


^1 

I  s 

Is 


5^  5  -Cj 

T**  ^  ^ 

»  c  v: 

^  C:  ^ 

5 

*>>  c: 

^  5  'C. 


Spanwise  station,  ^ 
(c)  Dead -weight  shear. 


Figure  5.-  Continued . 


Aeroelastic  deflection  per  unit  Aeroeiastic  loading  coefficient  per 

amber  (c'^^  ^-Oi),  e^,  ,  rad  unit  camber  (c^^^-.OI),  Cj^  c/Cg, 


Span  wise  station,  y 
(d)  Torsional  loading  due  to  camber. 
Figure  5.-  Concluded. 
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Dynamic  pressure,  ib/sq  ft 


Figure  6.-  Variation  in  lift-curve-slope  ratio  with  dynamic  pressure  for 
example  wing. 
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